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ABSTRACT. The first molecular phylogeny is presented for four species of the scorpion genus Meso- 
buthus , based on DNA sequences of three gene fragments (two mitochondrial and one nuclear protein 
coding gene, ~1 kb). The inferred phylogeny based on a pooled maximum likelihood analysis revealed a 
clear deep splitting between the “western clade” consisting of M. gibbosus and M. cyprius (Greece/ 
Anatolia, Cyprus) and the “eastern clade” consisting of M. eupeus and M. Caucasians (Anatolia/Central 
Asia). The species M. Caucasians (recently placed in the genus Olivierus Farzanpay 1987) groups mon- 
ophyletically within Mesobuthus ; thus, the genus Olivierus is synonymized here with Mesobuthus. Se¬ 
quences of M. eupeus and M. caucasicus sampled mainly from Kazakhstan and Uzbekistan are highly 
structured, indicating the possible existence of multiple species. 
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The genus Mesobuthus Vachon 1950 (Scor- 
piones, Buthidae) currently includes 12 spe¬ 
cies (Fet & Lowe 2000; Gantenbein et al. 
2000b); its type species is Mesobuthus eupeus 
(C.L. Koch 1839). Except for M. gibbosus 
(Brulle 1832), which is found in the Balkans 
and Turkey, the diversity of this genus is con¬ 
centrated in Asia. Numerous species and sub¬ 
species are distributed from Turkey to Korea, 
with the centers of diversity in Central Asia 
and Iran. Mesobuthus species are the most 
common and abundant scorpions in a variety 
of arid habitats, from sand deserts to high 
mountains over 3000 m (Fet 1989, 1994). 
They are found up to 50°N in Kazakhstan, the 
northern limit of the natural range of the Old 
World scorpions (Gromov 2001). 

Although Buthidae are the most diverse and 
medically important family of scorpions (Fet 
& Lowe 2000), there has been no attempt so 
far to produce a phylogenetic analysis of this 
family. Especially powerful are phytogenies 
based on DNA sequence data in combination 
with morphology (Gantenbein et al. 1999a, 
2000a; Fet et al. 2001). Our first applications 
of this technique in Buthidae refer to the gen¬ 


era Buthus Leach 1815 (Gantenbein et al. 
1999b; Gantenbein & Largiader 2003) and 
Centruroides Marx 1890 (Gantenbein et al. 
2001; Towler et al. 2001); we also published a 
pilot phylogeny of 17 buthid genera (Fet et al. 
2003). Molecular markers helped to define is¬ 
land species, where neither the biological spe¬ 
cies concept nor any other species concept can 
be applied (Gantenbein et al. 2000b, 2001). 

METHODS 

The currently studied available material be¬ 
longed to four species: Mesobuthus gibbosus 
(Greece, Turkey), M. cyprius Gantenbein & 
Kropf 2000 (Cyprus), M. eupeus (C.L. Koch 
1839) (Turkey, Kazakhstan, Turkmenistan, 
Uzbekistan, China), and M. caucasicus (Nord- 
mann 1840) (Kazakhstan, Turkmenistan, Uz¬ 
bekistan) (see Table 1 for locality informa¬ 
tion). For DNA analyses, the total DNA was 
extracted from fresh or preserved (94-98% 
ethanol) muscle tissue using a standard phe¬ 
nol/chloroform and precipitation method 
(Sambrook et al. 1989). We amplified a ca. 
450 base pair (bp) fragment of the 16S rRNA 
mitochondrial (mt) DNA by polymerase chain 
reaction (PCR) using the primers and condi- 
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Table 1.—Sampling sites and country of origin of Mesobuthus species used in this study. 


Abbreviation Country 

Locality 

EMBL accession nos. 

Coll. 16S coxl PK 

M. caucasicus (Nordmann 1840) 




McrzKZal Kazakhstan 

Kapchagai 

Ch. Tarabaev, AJ550674 
1990 

AJ550692 

AJ550713 

McaKZb 1 Kazakhstan 

Baigakum 

V. Fet & A. AJ550675 

Gromov, 

25-V-2002 

AJ550693 

AJ550714 

McaUZal Uzbekistan 

Bukhara 

V. Fet & A. AJ550676 
Gromov, 

2G-IV-2002 

AJ550694 

AJ550715 

McciUZb 1 Uzbekistan 

Jarkurgan 

V. Fet & A. AJ550677 

Gromov, 

26-IV-2002 

AJ550695 

AJ550716 

McaUZcl Uzbekistan 

Karakalpak 

Steppe, 

Ferghana 

V. Fet & A. AJ550678 

Gromov, 

20-V-2002 

AJ550696 

AJ550717 

McaXJZd ! Uzbekistan Babatag 

Mountains 

M. cyprius Gantenbein & Kropf 2000 

V. Fet & AJ550679 

A. Gromov, 

30-IV-2002 

AJ550697 

AJ550718 

McyC Yal Cyprus 

Tepebasi 

A. Scholl, 27- AJ550680 
IX-1997 

AJ550698 

AJ550719 

McyCYb I Cyprus 

M. eupeus mongolicus (Birula 

Kantara 

1911) 

A. Scholl, 20- AJ550681 
V-1998 

AJ550699 

AJ550720 

MeuGobil China 

M. eupeus eupeus (C.L. Koch 

Gobi Desert 

1839) 

A. Davidson, AJ550682 

30-VII- 
1998 

AJ550700 

AJ550721 

McnTRa 1 Turkey 

Giilsehir, 

Central 

Anatolia 

A. Scholl, 28- AJ550688 
V-1998 

AJ550701 

AJ550722 

MeuTRb 1 -2 Turkey 

Cemilkoy, 

A. Scholl, 31- AJ550689- 

AJ550702- 

AJ550723- 

Central 

Anatolia 

M. eupeus thersites (C.L. Koch 1839) 

V-1998 90 

03 

24 

MeuKZal Kazakhstan 

Bakanas 

A. Gromov, AY228141 

2-5-VI- 
2000 

AJ550704 

AJ550725 

MeuKZbl Kazakhstan 

Baigakum 

V. Fet & A. AJ550684 
Gromov, 

23-V-2002 

AJ550705 

AJ550726 

MeuKZcl Kazakhstan 

Karatau 

Mountains 

V. Fet & A. AJ550685 

Gromov, 

27-V-2002 

AJ550706 

AJ550727 

MeuTU 1 a Turkmenistan 

Repetek, 

Karakum 

V. Fet & A. AJ550686 
Gromov, 

15-18-IV- 

2002 

AJ550707 

AJ550728 

MeuUZ&l Uzbekistan 

Zarafshan, 

Kizylkum 

A. Gromov, AJ550687 

18-IV-1998 

AJ550708 

AJ550729 
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Table I.—Continued. 


EMBL accession nos. 


Abbreviation 

Country 

Locality 

Coll. 

16S 

cox I 

PK 

MeuU Zbl 

Uzbekistan 

Babatag 

V. Fet & A. 

AJ550683 

AJ550709 

AJ550730 



Mountains 

Gromov, 

29-IV-2002 




Mesobuthus gibbosus (Brulle 1832) 





Mg/GRal 

Greece 

Mathia, Pelo- 

I. & B. Gan¬ 

AJ402571 

AJ550710 

AJ550731 



ponnesos 

tenbein, 18- 
III-1998 




Mg/GRb 1 

Greece 

Igoumenitsa, 

I. & B. Gan¬ 

AJ550691 

AJ550711 

AJ550732 



Epirus 

tenbein, 28- 
IV-1998 




Mg/TRal 

Turkey 

Avanos, Cen¬ 

A. Scholl, 28- 

AJ402587 

AJ550712 

AJ550733 



tral Anato¬ 
lia 

V-1998 




Outgroup: Androctonus australis (Linnaeus 1758) 




AanTNal 

Tunisia 

Nefta 

A. Scholl, 27- 

AJ506868 

AJ506919 

AJ550734 




IV-1999 





tions of Gantenbein et al. (1999). For the par¬ 
tial amplification of the cytochrome oxidase I 
(carl) gene we used the primers LCO (Folmer 
et al. 1994) and Nancy (Simon et al. 1994) 
which amplify a —850 bp fragment (positions 
—30-850 of the Drosophila yakuba sequence; 
Flybase: FBgn0013179). We used the follow¬ 
ing PCR profile: initial denaturation 94 °C for 
4 min followed by 40 cycles of 25 s at 94 °C, 
20 s at 51 °C, and 90 s at 72 °C. In addition, 
we designed new buthid-specific PCR primers 
from a clone (03B09) of an EST library of 
Mesobuthus gibbosus (unpublished data, Gan¬ 
tenbein et al., in preparation) to amplify a 
—360 bp fragment of the nuclear protein ki¬ 
nase (PK) gene (Flybase locus CGI 1221, 
identified from BLASTX against the Dro¬ 
sophila protein database, similarity 43%, Ex¬ 
pectation = 5e-14). The primers were 
03B09for 5'-TCT GAT GTA TGG CAG ATG 
GCA ATG-3' and 03B09rev 5'-CGA ACT 


CAA GAT CCA CTC CTG TAC TCG-3'. We 
used the same PCR profile as for carl. PCR 
primers were removed by polyethylene glycol 
(PEG 8000) precipitation, and templates were 
directly sequenced on one strand using one of 
the PCR primers and DYEnamic ET Dye Ter¬ 
minator Kit (Amersham Biosciences). For 
carl and PK we used the forward PCR-prim- 
ers for sequencing, and for the 16S we used 
the same primer as in Gantenbein et al. 
(1999). Sequencing reactions were ethanol/so¬ 
dium-acetate precipitated and run on an 
ABI377XL sequencer (Applied Biosystems, 
Foster City, CA). All sequences were checked 
manually for sequencing errors. All sequences 
were deposited in the EMBL nucleotide se¬ 
quence database (http://www.ebi.ac.uk). As 
outgroup, we used the Old World buthid An - 
droctonus australis (L. 1758) from Nefta, Tu¬ 
nisia (AauTNal), which is a suitable outgroup 
for the genus Mesobuthus as indicated from 


—> 

Figure 1.—Maximum Likelihood (ML) tree of Mesobuthus species from southern Europe, Western and 
Central Asia inferred from three combined DNA sequence fragments of the mitochondrial 16S, cox I and 
the nuclear PK regions, 1,095 bp (-In Likelihood was 5505.86). The DNA substitution model was TRN 
+ T + I (Tamura & Nei 1993); base frequencies: i r A = 0.29, tt t = 0.37, ir c = 0.14, R matrix = (A G = 
6.70, A-C = A-T = G-T = 1, C-T = 3.56), gamma shape parameter a = 0.53, and proportion of invariable 
sites = 0.37, respectively. The tree was rooted using the outgroup species Androctonus australis (Aau¬ 
TNal). Numbers at nodes refer to bootstrap support given as percentage from 1,000 pseudo replicates by 
neighbour-joining of ML distances. 
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Table 2.—Estimates of average rates of synonymous substitutions ( K s ) (lower left) and non-synonymous 
(K a ) (upper right) substitution per site among four species of Mesobuthus of two protein-coding gene 
fragments. Diagonals show average K s and K A (in parentheses) among haplotypes within species. Estimates 
of nucleotide diversity rr for each species (Nei & Li 1979) are given in the last column. 



cyprius 

caucasicus 

eupeus 

gibbosus 

diversity it 

coxl locus 

cyprius 

0.11 (0.009) 

0.297 

0.024 

0.006 

0.037 

caucasicus 

0.844 

0.482 (0.030) 

0.034 

0.023 

0.098 

eupeus 

0.887 

0.753 

0.377 (0.026) 

0.023 

0.089 

gibbosus 

0.538 

0.892 

0.692 

0.223 (0.004) 

0.048 

PK locus 

cy prius 

0.038 (0.018) 

0.106 

0.013 

0.013 

0.032 

caucasicus 

0.093 

0.022 (0.000) 

0.004 

0.004 

0.008 

eupeus 

0.073 

0.032 

0.008 (0.000) 

0.009 

0.002 

gibbosus 

0.04 

0.082 

0.064 

0.037 (0.009) 

0.017 


morphology (Vachon 1952). Voucher speci¬ 
mens will be deposited in the Natural History 
Museum, Bern, Switzerland (NHMBE). 

The resulting 22 16S DNA sequences were 
aligned using ClustalX (Thompson et al. 
1997) and by eye. We applied Maximum 
Likelihood (ML) to the DNA sequence data. 
The alignment was 344 bp long. The align¬ 
ment was deposited in the EMBL nucleotide 
database (ALIGN_000522). Ambiguities and 
gaps were stripped out (Swofford et al. 1996), 
leaving 302 bp. The alignment of the coxl and 
the PK fragments was unambiguous because 
of the open reading frames (ORF). In order to 
choose the most appropriate DNA model of 
nucleotide substitution, we calculated hierar¬ 
chic likelihood ratio test statistics using MO- 
DELTEST 3.06 (Posada & Crandall 1998; 
Huelsenbeck & Rannala 1997) which is im¬ 
plemented in PAUP* 4.0b 10 (Swofford 1998) 
and calculates the hierarchical likelihood ratio 
statistics (LRT) of 56 different substitution 
models based on a NJ tree using JC69 dis¬ 
tances (Jukes & Cantor 1969). The rate het¬ 
erogeneity among sites was assumed to follow 
a gamma distribution (shape parameter a was 
ML-estimated) with four categories, each rep¬ 
resented by its mean (Yang 1996). The ML- 
estimated parameters and the model are given 
in the legend of Fig. 1 

The tree topology found with ML of the 
pooled mitochondrial data (16S & cox I) was 
not significantly different from the topology 
inferred from the nuclear PK gene if assessed 
by the two-tailed K-H-test (-lnL (treePK | PKdata) = 
618.49, lft^(treel6S&cojd|Pkdata) 632.50, P 


0.16) with re-estimation of maximum likeli¬ 
hood by non-parametric bootstrapping 
(RELL) (Kishino & Hasegawa 1989; Kishino 
et al. 1990). Thus, we pooled the nuclear and 
mitochondrial data, leaving 1,095 bp (475 bp 
of the cox I and 318 bp of the PK, excluding 
gaps and ambiguities). 

In a further step, the molecular clock hy¬ 
pothesis (i.e., equal rates across all sequences) 
was tested using the x 2 approximated likeli¬ 
hood ratio test statistics with OTU’s—2 de¬ 
grees of freedom (df = 22-2 = 20) which was 
rejected with a P -value < 0.01 (Huelsenbeck 
& Crandall 1997). We explored the tree space 
by 100 heuristic tree searches using the Tree- 
Bisection-Reconnection (TBR) algorithm and 
by randomizing the order of the sequence in¬ 
put in PAUP*. Phylogenetic trees were rooted 
using Anclroctonus australis as an outgroup. 
Statistical confidence of phytogenies was as¬ 
sessed using the bootstrap procedure (1,000 
pseudoreplicates) (Felsenstein 1985) using 
PAUP*. These distances were usually —15 % 
between Mesobuthus species and go up to 
—20 % if Mesobuthus sequences were com¬ 
pared to the outgroup. 

RESULTS AND DISCUSSION 

The three analyzed DNA fragments con¬ 
tained considerable polymorphism among but 
also within species as estimated using Nei & 
Li’s (1979) nucleotide diversity i r and the 
amount of synonymous versus non-synony¬ 
mous substitutions of coding regions (Jukes & 
Cantor 1967) (Table 2). We found 153 muta¬ 
tions out of 475 bp (of which 29 were replace- 
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ment changes) in the cox I fragment, and 19 
mutations out of 318 bp (of which 8 were re¬ 
placement changes) in the PK gene (excluding 
the outgroup). The 16S fragment contained 
131 polymorphic sites out of 302 bp. 

The recovered phylogeny (Fig. 1) showed 
high support for all four included species. 
This phylogeny allows to address several im¬ 
portant taxonomic and evolutionary issues 
pertaining to the genus Mesohuthus. First of 
all, the phylogeny demonstrates a deep split 
between the “western clade” of M. gibbosus 
and M. cyprius and the “eastern clade” of M. 
eupeus and M. caucasicus (Fig. 1); each clade 
is well supported (100% and 92% bootstrap, 
respectively). Within the “eastern clade,” 
there is a strong support for currently accepted 
species M. eupeus (87%) and M. caucasicus 
(94%). The Mesobuthus caucasicus sequences 
are nested within the genus Mesobuthus as a 
sister group to M. eupeus. This observation is 
important as the generic placement of M. cau¬ 
casicus has been controversial since Farzan- 
pay (1987) created a separate monotypic ge¬ 
nus Oiivierus for this species. As our 
molecular data show, this new genus is para- 
phyletic with respect to Mesobuthus (Fig. 1). 
Fet & Lowe (2000) listed Oiivierus as a valid 
genus but noticed that it was created without 
any solid justification or revision. The only 
characters (number of granules on movable 
finger of pedipalp chela) that Farzanpay 
(1987) briefly quoted (in Farsi) in support of 
Oiivierus , were borrowed from Vachon (1958) 
species-level descriptions, and are not diag¬ 
nostic for a genus. Further, our molecular data 
do not support monophyly of Oiivierus. 
Therefore, we propose to list this genus as a 
synonym of Mesobuthus , as was traditionally 
accepted by all authors before Farzanpay 
(1987): Oiivierus Farzanpay 1987 = Meso¬ 
buthus Vachon 1950, NEW SYNONYMY. 

Another observation refers to the subspe¬ 
cific structure of M. eupeus , the type species 
of the genus Mesobuthus. This species is ex¬ 
tremely polymorphic, and has 14 formally val¬ 
id subspecies ranging from Turkey to China 
(Vachon 1958; Fet 1989; Fet & Lowe 2000), 
with most subspecies described from Iran. The 
nominotypic subspecies M. e. eupeus (C.L. 
Koch 1839) is found in the Caucasus and Tur¬ 
key, while most populations from Central 
Asia are classified as M. e. thersites (C.L. 
Koch, 1839). In our analysis, the nominotypic 


subspecies (Turkey) is highly supported (boot¬ 
strap 100%) while the Central Asian sequenc¬ 
es present a more complicated case. A sepa¬ 
rate clade exists for a sand desert population 
from Turkmenistan (MeuTUal), as opposed to 
the clade of several populations from Uzbek¬ 
istan, Kazakhstan and China (bootstrap 68%). 
This can be an indication of an ancient sepa¬ 
ration between southern and northern desert 
forms, possibly valid at least at the subspecies 
level (Fet 1994). The further phylogenetic 
analysis could result in elevating these sub¬ 
species to the species level. It is interesting 
that already Birula (1917) grouped all subspe¬ 
cies of M. eupeus into two species groups 
(“sections”), “ eupeus ” and “thersites”; how¬ 
ever, status of these groups was never exam¬ 
ined. 

Within M. caucasicus , a very strongly sup¬ 
ported clade (bootstrap 98%) groups popula¬ 
tions from central Uzbekistan (Bukhara, Fer¬ 
ghana) and Kazakhstan, while those from 
southern Uzbekistan (Jarkurgan, Babatag) 
group outside. It remains to be seen if genetic 
separation in this case is matched by the mor¬ 
phological variation, as there are several sub¬ 
species described from Central Asia as well. 

In addition, in our phylogeny the northern 
Central Asian populations of both Mesobuthus 
eupeus and M. caucasicus (Kazakhstan) ap¬ 
pear to be derived compared to the southern 
populations of both species (Uzbekistan); this 
could be the result of progressive Tertiary ar- 
idization and spreading of the arid scorpion 
species from south to north to the sand and 
clay deserts (Fet 1994). 

The presented data also allow a calibration 
of a molecular clock using the separation of 
Cyprus from the Anatolian mainland (5.2 
Mya) after the Messinian salinity crisis, dur¬ 
ing which gene flow between island and main¬ 
land populations could have been possible. 
The Mediterranean Basin was refilled within 
only 100 yrs, which provides an excellent cal¬ 
ibration point for a molecular clock (Hsu et 
al. 1977; Gantenbein & Largiader 2002). 
Thus, the sequence divergence between Ana¬ 
tolia and Cyprus was estimated to 0.09 ± 0.01 
(0.10 ± 0.01 for 16S), which results in a se¬ 
quence divergence rate of 0.017 per My. This 
rate estimate is somewhat higher than previ¬ 
ous estimates in scorpions for Mesobuthus 
gibbosus (Gantenbein & Largiader 2002) but 
lies absolutely in the range of scorpions such 
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as Buthus occitanus and Centruroides (Buth- 
idae) and other invertebrates such as butter¬ 
flies, beetles and crickets (Brower 1994; 
Fleischer et al. 1998; Gomez-Zurita et aL 
2000; Gantenbein et ah 2001; Gantenbein & 
Largiader 2003). 

The genus Mesobuthus was created when 
Vachon (1950) initiated a large-scale “split¬ 
ting” revision of the traditional genus Buthus 
Leach 1815. Its composition is still contro¬ 
versial. For instance, several Indian species 
are classified currently in this Mesobuthus but 
their generic identity is unclear (Fet & Lowe 
2000). Separate species have been only re¬ 
cently described from Cyprus (Gantenbein et 
aL 2000b) and confirmed for northern Israel 
(Fet et al. 2000), Numerous subspecies exist 
in M. eupeus and M. caucasicus but morpho¬ 
logical characters are inconclusive as for their 
relationship and taxonomic status. Further ap¬ 
plication of new molecular markers will fa¬ 
cilitate our understanding of taxonomy and 
evolution of this common scorpion genus. 
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